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We employ Optimal Control Theory to discover an efficient information retrieval algorithm that 
can be performed on a Rydberg atom data register using a shaped terahertz pulse. The register is a 
Rydberg wave packet with one consituent orbital phase-reversed from the others (the "marked bit" ) . 
The terahertz pulse that performs the decoding algorithm does so by by driving electron probability 
density into the marked orbital. Its shape is calculated by modifying the target of an optimal control 
problem so that it represents the direct product of all correct solutions to the algorithm. 
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I. INTRODUCTION 

The interaction of terahertz frequency electromagnetic pulses with Rydberg atoms has produced many insights 
into the dynamical properties of atomic systems [Q. The comparable time scales of terahertz pulses with those of 
Rydberg state lifetimes make it possible to envision schemes of quantum control. In this paper, we propose a method 
of controlling Rydberg wave packets using shaped terahertz pulses, and theoretically show how these pulses can be 
designed to execute a quantum algorithm on a Rydberg atom data register. 

It has been shown that information can be stored in the phases of the constituent orbitals of a Rydberg wave 
packet [HrI . Recently, a terahertz half-cycle pulse was used to decode the information stored in a Rydberg atom data 
register 0^ This half-cycle pulse decodes the phase structure by retaining the population only in the orbital that was 
initially 180° out of phase with respect to the other orbitals, i. e. , the marked bit. However, this guess (unshaped) 
pulse does not decode all marked bits of the register with the same efficiency. In this paper, we aim to find the shaped 
THz pulse that will optimally transfer most of the population to any marked bit of the quantum data register. 

II. OPTIMAL CONTROL THEORY 

To design the terahertz pulse, we use a method that has been used extensively in mathematical and engineering 
applications — Optimal Control Theory (OCT) 1^,^. This theory has also been applied to the control of quantum 
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systems [|6|^2| with some success in experimental implementation 
pulse that can be used to achieve a desired target state. We then modi 
discover not a single target, but an optimized quantum algorithm. The wave function of the Rydberg electron is the 
state variable, and the electric field of the terahertz pulse is the control parameter. A functional J is defined, whose 
extremum must be calculated. The functional consists of two parts, representing the desired target and the cost. 
Our aim is to maximize the fraction of the electron probability density in a target orbital |a/j) at a time T (after the 
end of the terahertz pulse). That is, the target functional {Pk{T)) — {ip{T)\ak) {ak\'ip{T)) must be a maximum. The 
cost functional represents the constraint on the control parameter, the terahertz field E{t). The integrated energy of 
the pulse must be kept low, therefore the cost functional is defined as Y{T) = dt £{t)\E{t)\'^ . Here ^ is a penalty 
parameter, in general time dependent, that controls the cost functional, and hence the peak terahertz field. The 
functional J written as 



j = (V'(r)|p,|v(T)}-r(r), (i) 

must be maximized. 

The aim is to find an optimal control function E{t) that maximizes J. The wave packet evolution is governed by 
the Schrodinger equation. In atomic units, e = rrig = h = 1, 

\m) = -iH{t)m)), (2) 

where H{t) = Hq + E{t)z. The equation of motion acts as a constraint on the evolution of the state, and in a 
manner similar to that used in variational calculus, we introduce a Lagrange multiplier \\{t)) p^ . The unconstrained 
functional that must be optimized is written as 
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Several iterative techniques for determining the optimal solution have been developed ]15|- |l7|jl9| , |lq | . Following the 
scheme for an iterative solution proposed in Ref. |p9[|, the functional J is written as the sum of a terminal part and 
an integral 

J = G+ [ dt R, where (5) 
Jo 

G - {^{T)\Pk\m)) - 2Re{Xit)m))\I , (6) 

R = -e{t)\E{t)\'' + 2Re [(A(t)lV(t)) - i{m\H\m)] ■ (7) 

The maximum of both G and R is sufhcient to ensure the maximum of J. Our objective is to iteratively determine 
the optimal function E{t) that maximizes J. The functional J at the fc*'* and (fc + 1)** iteration is defined by the 
changes in the wave function and the field. That is, 
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The difference in J between two successive iterations is written as 

jfc+i _ Jfe ^ Ai + A2 + A3, where (9) 
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|A(t)) = -ii?'=|A(t)},we find (14) 
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A3 = 0. (17) 

The solution improves when J increases or stays the same at each iteration. Pk is a positive semidefinite operator, 
therefore Ai is greater than or equal to zero. The change in the control parameter AE is chosen to maximize A2. 
The expression for A2 suggests the following change in the field at the (k+l)th iteration p9| ]: 

AE{t)^j^^{X'^{t)\z\i,'^+\t)}. (18) 

The appearance of 1^^^+^ (t)) in the above expression implies that the overlap of \ip(t)) and \ X{t)) is fed back immediately 
to find the field at the next time step. 

The optimal control algorithm consists of the following steps: 

1. Starting from the initial wave packet |?/'^°-'(0)) = |'0(O)) and a first guess for the terahertz field E'-^\t), the wave 
packet is propagated according to Eq. || to find \^/^'>{T)). 

2. Using Eq. |l^ to find |A(r)), Eq. |l^ is iterated backward to time t — 0, and |A(t)) is found at every time step. 

3. With |'i/;^(0)) = \ip{0)), equation ^s] is then used to find a new value of the control field E^{t) and equation || is 
used to propagate the wavepacket forward in time. 

The second two steps are repeated until the target yield converges to within the desired accuracy. 



III. OPTIMAL PULSE FOR A SINGLE TARGET STATE 



The first step in our approach is to find the optimal THz field required to decode a single flipped state. The initial 
state of the Rydberg data register is a wavepacket made of the 2Ap through 29p orbitals of equal amplitudes and the 
phase of the 26p (the marked bit) orbital opposite to that of the others. The initial guess THz pulse is a half-cycle 
pulse of pulse width Ips. The desired initial phase structure occurs at the peak of the half-cycle pulse (i. e. at 
0.5ps). We find the THz pulse that will optimally transfer most of the population to the marked bit. The best value 
of the penalty parameter, i, that controls the peak field of the THz pulse to a reasonable value, and at the same 
time produces the desired final state, was found to be roughly lO^*'. One feature of several optimal fields obtained 
theoretically |0,|ll| is that the fields do not go to zero smoothly at the times t = and t = T. In the present system, 
the correct evolution of the wave packet depends very sensitively on the fields at the end points. The addition of a 
smooth switch-on and switch-off of the calculated optimal field drastically changes the evolution of the wave packet. 
Therefore, the condition that the THz field continuously goes to zero before and after the time interval of choice must 
be built into the algorithm. To ensure that the THz field goes smoothly to zero at times t = Q and t — T, the penalty 
parameter £[t) is made a smoothly varying time-dependent function. The penalty on the pulse fluence is a thousand 
times more at the end points than at the rest of the pulse duration. The smoothness of the penalty function ensures 
a smooth switching on and switching off of the THz pulse. 

This OCT implementation is very successful in describing the terahertz control of a Rydberg wave packet. The 
method takes macrosteps in the control field at every iteration and convergence is swift. The computational complexity 
is of the same order as that of the wave packet propagation. Therefore, we use a split-operator method in a restricted 
basis of essential stat es pl| ,^ . The energy eigenstates of cesium are calculated using a pseudopotential method on a 
nonlinear radial grid ]22|]. A time step of lOfs ensures the accuracy of the propagator, which is correct through the 
second-order in the time step. The numerical implementation of the local iterative algorithm is extremely sensitive 
to numerical error, and A3 must be maintained equal to zero to very high precision ||T^. The unitary nature of the 
symmetrized-product propagator maintains this condition. The restricted basis consists of 195 energy eigenstates 
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with principal quantum number n ranging between 21 and 31, and angular momentum quantum number £ < 17. 
Absorbers ensure that population does not get reflected from the n — 21, n = 31, and £ = 16 'boundaries'. Using the 
selected state basis also imposes the condition that the spectrum of the THz pulse is decided by the energy range of 
the selected state basis. In this implementation of optimal control theory, we have chosen a fixed pulse length, T, of 
roughly 8ps. In other formalisms, this time T may also be varied as a parameter. 

The THz field that optimizes the population in the marked 26p state is shown in Fig. ??(a). The initial population 
in the 26p state is 16.7%. The optimal pulse will decode the information stored as phase by transferring most of the 
population into the 26p state. With the initial guess pulse, the population is 29.5%. After 50 iterations, the target 
yield is increased to 52.8%. The spectrum and Husimi distribution Q of this optimal pulse are shown in Fig. ??(b) 
and Fig. ??(c) respectively. Notably, the strong peaks in the spectrum and the Husimi distribution do not correspond 
to any resonance between the energy levels of the selected state basis. The optimal terahertz pulse does not drive 
the system to any particular resonant condition. Instead, it alters the phases of the constituent orbitals of the wave 
packet so that they interfere to produce the desired probability distribution. 

Figure ?? shows the evolution of the wavepacket as a function of time while the optimal pulse is on. During the 
THz pulse, probability can leak into other states not in the register (the other states in the essential basis). At the 
end of the pulse, a large fraction of the electron probability density lies in the flipped orbital (marked bit) of the data 
register. This can be thought of as using the other states of the data register as working qubits, which are used during 
the computation, but are not measured for any useful retrieval of information. 

One interesting feature of this optimal pulse is that the peak field of roughly IKV/cm lasts for roughly 0.5ps. For 
a n = 26 wavepacket, this field which is beyond the field ionization limit lasts for more than half the Kepler period 
(- 27rn3). Yet, 99% of the population remains in the selected state basis. This feature is an example of interferometric 
stabilization p^], seen in other atomic systems. 



IV. OPTIMAL PULSE FOR A QUANTUM ALGORITHM 



This THz pulse is optimal only for decoding the flipped 26p orbital. That is, if the phase of a different state were 
flipped, this pulse will not decode it. We wish to design a universal THz pulse that will optimally decode any flipped 
orbital of the wave packet register. Therefore we redefine the optimal control problem by considering an initial state 
that is a product state of independent wave packets with singly flipped orbitals. 

i*(o)) = Hglm ® V'feiw) ® 4%m ® • • • (19) 

The terahertz pulse acts simultaneously, but independently on all these wave packets. The desired final state is also 
a product state of independent wave packets with the fiipped bit correctly decoded. 

|4'(r)) = |25p(i)) ® |26p(2)) |27p(3)) ® • • • (20) 

The counterparts of Eq. ^ and Eq. ^ are straightforward. At every time step, the updated THz field is found by 
using a modified version of Eq. ^ with the matrix element of z replaced by a sum of matrix elements of z, one from 
each independent 'subspace'. 

-1 ^ 

= ^E(4)Wi*^i^w'w)- (21) 

Using this method, we find the terahertz pulse that detects any flipped orbital of the A^-bit data register. The 
advantage of this reflnement is that the computational resources needed increase only by a factor of the number of 
constituent states in the wave packet register. 

We now flnd the optimal terahertz pulse that will decode any flipped state in a six state Rydberg data register. 
The register consists of np states of cesium, with n from 24 — 29. Population in the flipped orbital is amplified by the 
diffusion of probability density from the adjacent states. This is an example of the implementation of Grover's search 
algorithm, where information is stored in states with differing phases, and a marked bit is amplified by "quantum 
diffusion" |25j. The outer states n — 2A and n = 29 are therefore not included in the optimization. The universal 
decoding pulse and its effect on a wave packet with different marked bits is shown in Fig. ??. After the pulse, the 
wave packet population is distributed so that the flipped state is clearly amplified. 
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V. CONCLUSIONS 



In conclusion, we have designed a terahertz pulse to implement a search algorithm on a quantum data register. 
Phase information stored in a Rydberg wave packet was optimally retrieved through the interaction with the pulse. 
Careful attention was paid to the smooth switch on and switch off of the THz pulse. We also show that it is possible 
to design an optimal pulse that can achieve not only a desired target state of an atom, but also implement a desired 
algorithm. 

To our knowledge, this is the first time that optimal control theory has been applied to the terahertz control of a 
quantum system. This theoretical study motivates the experimental design and control of terahertz frequency pulses. 
Beyond quantum control, these results point to the possibilities of using Rydberg atoms as quantum computers, and 
terahertz pulses to implement quantum algorithms. 
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